set more off

use dataset, clear

local Z pcteligible2013

keep if abs(x) < 100

su `Z', det
gen Z = `Z' >=r(p50)
replace Z = . if mi(`Z')

su demvoteshare, det
local mymed = r(p50)

preserve

	keep if demvoteshare < `mymed'

	#delimit

	gr tw 
		(sc dturnout20142010 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dturnout20142010 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20142010 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g1, replace)
			;		
		
	#delimit cr

	#delimit

	gr tw 
		(sc dturnout20142010 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dturnout20142010 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20142010 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("High eligibility")
			legend(off)
			nodraw
			name(g2, replace)
			;
		
	#delimit cr


	#delimit

	gr tw 
		(sc dturnout20162012 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dturnout20162012 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20162012 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g3, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dturnout20162012 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dturnout20162012 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20162012 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("High eligibility")
			legend(off)
			nodraw
			name(g4, replace)
			;
		
	#delimit cr

	gr combine g1 g2, rows(1) title("Change in turnout 2014-2010") name(G1, replace) ycommon
	gr combine g3 g4, rows(1) title("Change in turnout 2016-2012") name(G2, replace) ycommon
	gr combine G1 G2, rows(2) title("Republican counties")

	gr export "_output/figureA18.pdf", replace

restore

preserve

	keep if demvoteshare > `mymed'

	#delimit

	gr tw 
		(sc dturnout20142010 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dturnout20142010 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20142010 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g1, replace)
			;		
		
	#delimit cr

	#delimit

	gr tw 
		(sc dturnout20142010 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dturnout20142010 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20142010 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("High eligibility")
			legend(off)
			nodraw
			name(g2, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dturnout20162012 x if Z == 0, msize(tiny) col(gs13))
		(lpolyci dturnout20162012 x if x < 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20162012 x if x > 0 & Z == 0, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("Low eligibility")
			legend(off)
			nodraw
			name(g3, replace)
			;
		
	#delimit cr

	#delimit

	gr tw 
		(sc dturnout20162012 x if Z == 1, msize(tiny) col(gs13))
		(lpolyci dturnout20162012 x if x < 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		(lpolyci dturnout20162012 x if x > 0 & Z == 1, clcol(black) lwid(medthick) ciplot(rline) alpat(solid) alcol(black) alwid(thin))
		,
			ylab(-20(10)10 , angle(horiz) nogrid)
			plotregion(style(none))
			subtitle(, fcol(none) lwid(none))
			xline(0)
			ylab(, angle(horiz))
			xtitle("Distance from nearest state border")
			legend(off)
			ytitle("Change in voter turnout")
			title("High eligibility")
			legend(off)
			nodraw
			name(g4, replace)
			;
		
	#delimit cr

	gr combine g1 g2, rows(1) title("Change in turnout 2014-2010") name(G1, replace) ycommon
	gr combine g3 g4, rows(1) title("Change in turnout 2016-2012") name(G2, replace) ycommon
	gr combine G1 G2, rows(2) title("Democratic counties")

	gr export "_output/figureA19.pdf", replace

restore
gr close
